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Abstract 

Data on solar and atmospheric neutrinos provide evidence for neutrino masses 
and mixings. We review some basic neutrino properties, the status of neutrino 
oscillations and implications for grand unified theories, including leptogenesis. 
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1 Neutrinos and Symmetries 



Neutrinos have always played a special role in physics due to their close connection with 
fundamental symmetries and conservation laws. Historically, this development started 
with Bohr who argued in 1930 that 

• energy and momentum 

might only be statistically conserved [T] in order to explain the continuous energy 
spectrum of electrons in nuclear /3-decay. As a 'desparate way out', Pauli then 
suggested the existence of a new neutral particle @, the neutrino, as carrier of 
'missing energy'. The emission of neutrinos in /3-decays could also explain the 
continuous energy spectrum. 

• lepton number 

As a massive neutral particle, the neutrino can be equal to its antiparticle || and 
thereby violate lepton number, a possibility with far reaching consequences. Most 
theorists are convinced that neutrinos are Majorana particles although at present 
there are no experimental hints supporting this belief. 

• parity and charge conjugation 

In the theory of weak interactions |4j chiral transformations [[| have played a crucial 
role in identifying the correct Lorentz structure. The resulting V — A theory H 
violates the discrete symmetries charge conjugation (C) and parity (P), but it 
conserves the joint transformation CP. 

• CP invariance 

The modern theory of weak interactions is the standard model with right- 
handed neutrinos. In its most general form CP invariance and lepton number are 
violated which, as we shall see, has important implications for the cosmological 
matter-antimatter asymmetry. 

• CPT and Lorentz invariance 

At present neutrino oscillations are discussed as a tool to probe CPT and Lorentz 
invariance || whose violation is suggested by some modifications of the standard 
model at very short distances. 

It therefore appears that we are almost back to the ideas of Bohr, with the hope for 
further surprises ahead. 



3 



2 Some Neutrino Properties 



Neutrino physics is a broad field || involving nuclear physics, particle physics, astro- 
physics and cosmology. In the following we shall only be able to mention some of the 
most important neutrino properties. Most of the lectures will then be devoted to neu- 
trino oscillations, the connection with grand unified theories and also the cosmological 
matter-antimatter asymmetry. 



2.1 Weak interactions 

A major source of neutrinos and antineutrinos of electron type is nuclear /3-decay, 



e > 



A(Z,N) -> A(Z + l,N-l) + e~ + u e 
A(Z,N) -> A(Z-l,N+l) + e + + v e , (1) 

which includes in particular neutron decay, n — > p+e~ +V e . Muon neutrinos are produced 
in 7r decays, 7r ± — > y£ + v fl (T' fl ), e ± + u e (V e ) and muon decays, — > e ± + V fl (u fl ) + u e (u e ). 
All these processes are described by Fermi's theory of weak interactions, 

L F = -^=J™J C ™ , (2) 

where Gf — 1.166 x 10 -5 GeV~ 2 QTO | is Fermi's constant. The charged current J^ c has 
a hadronic and a leptonic part, 

joa = jcc (h) + jcc (l) ( (3) 

J» C{h) = PJ»(9v-gAK)n + Ud tl 7r + + ... , (4) 
tf C{l) = ^7m(!- 75)6 + ^(1-75)// + ... . (5) 

Here /„■ is the pion decay constant, and gy = 0.98 and ^ = 1.22 are the vector and 
axial-vector couplings of the nucleon. 

An impressive direct test of the V — A structure of the charged weak current is the 
helicity suppression in tt decays. Since pions have spin zero, angular momentum con- 
servation requires the same helicities for the outgoing antineutrino and charged lepton. 
On the other hand, the V — A current couples particles and antiparticles of opposite 
helicities. Hence, one helicity flip is needed, and the decay amplitude is proportional to 
the mass of the charged lepton. This immediately yields for the ratio of the partial decay 
widths, 

r(7r -> fi + zaj v^w \ m i- m iJ 
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Figure 1: Kurie plot for tritium beta- decay. 



in remarkable agreement with experiment 



R exp = (1.230 ± 0.004) x 10" 



(7) 



From the measurement of the invisible width of the Z-boson we know that there are 
three light neutrinos [[LTJ, 

r 



N,, 



r,. 



2.984 ±0.008 . 



This is consistent with the existence of u e , and u T , i.e. one neutrino for each quark- 
lepton generation. Direct evidence for the r-flavour of the third neutrino has recently 
been obtained by the DONUT collaboration at Fermilab, which observed r-appearance 



in nuclear emulsions 11 



v T + N -> r + X . 



(9) 



The same reaction will be used to prove v„ — > v T oscillations by r-appearance. 



2.2 Bounds on neutrino masses 

Neutrinos are expected to have mass, like all other leptons and quarks. Direct kinematic 
limits for tau- and muon-neutrinos have been obtained from the decays of r-leptons and 
7r- mesons, respectively. The present upper bounds are [ jH| , 

m Ur < 18.2 MeV (95% CL) , < 170 keV (90% CL) . (10) 
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Figure 2: Effective electron mass in Qv2f3-decay. 

These are bounds on combinations of neutrino masses, which depend on the neutrino 
mixing. The study of the electron energy spectrum in tritium /3-decay over many years 
has led to an impressive bound for the electron-neutrino mass. The strongest upper 



bound has been obtained by the Mainz collaboration [12| 

m Ue < 2.2 eV (95% CL) . (11) 

It is based on the analysis of the Kurie plot (fig. |l|) where the electron energy spectrum 
is studied near the maximal energy Eq, 



K(E e ) oc J(E - E e )((E - E e y - m*)W 



(12) 



In the future the bound (|TT|) is expected to be improved to 0.3 eV [[13 



The most stringent bound so far has been obtained for the effective electron-neutrino 
mass m ee in neutrinoless double /3-decay (fig. 0), 



A(Z, N) -> A(Z + 2, N - 2) + e~ + e" . 



The Heidelberg-Moscow collaboration presently quotes the bound [[14 

< 0.38 k eV (95% CL) , 



(13) 



(14) 



where k = 0(1) represents the effect of the hadronic matrix element. The quantity m ee 
is a Majorana mass. Hence, a non-zero result would mean that lepton number violation 
has been discovered! The GENIUS project has the ambitious goal to improve the present 
bound (O) by about two orders of magnitude (15 . 



2.3 Cosmology and astrophysics 

Big bang cosmology successfully describes the primordial abundances of the light ele- 
ments 4 He, D and 7 Li. This also leads to predictions for the baryon density Qb and 
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Figure 3: Z-burst in the annihilation of ultrahigh energy cosmic neutrinos and relic neu- 
trinos. 

the effective number of light neutrinos at the time of neutrino decoupling (T ~ 1 MeV, 
t ~ 1 s) 

1.2 < N e J f < 3.3 (99% CL) . (15) 

Before the precise measurement of the Z-boson width this bound was already a strong 
indication for the existence of not more than three weakly interacting 'active' neutrinos. 
The bound also strongly constrains the possibility of large lepton asymmetries at the 
time of nucleosynthesis and the possible existence of 'sterile' neutrinos which do not 
interact weakly. Detailed studies show that mixing effects would be sufficient to bring a 
sterile neutrino into thermal equilibrium. This would imply N^** = 4, in contradiction 



to the bound fllS]) . The existence of a sterile neutrino is therefore disfavoured |T6 |. 

The standard cosmological model further predicts the existence of relic neutrinos with 
temperature T ~ 1.9 K and an average number density around 100/cm 3 per neutrino 
species. For neutrino masses between 3 x 1CT 2 eV (see section 4.3) and 2 eV (cf. ([TT|)) 
this yields the contribution to the cosmological energy density, 

0.001 < Vt u h 2 < 0.1 ; (16) 

here Q u = p u /p c , where p v is the neutrino energy density and p c = 3Hq/(37iGn) = 
1.05/i 2 x 10 4 eV/cm 3 is the critical energy density. Hence, Q u is substantially smaller 
than the contributions from the 'cosmological constant' and from cold dark matter 
(flcDAi), but it may still be as large as the contribution from baryons, VLsh 2 ~ 0.02. 

The detection of the relic neutrinos is an outstanding problem. It has been suggested 
that this may be possible by means of 'Z-bursts' (fig. |3|), which result from the resonant 
annihilation of ultrahigh energy cosmic neutrinos (UHEC^s) with relic neutrinos |18|, [19]. 
On resonance, the annihilation cross section is enhanced by several orders of magnitude. 
In the rest system of the massive relic neutrinos the resonance condition is fullfilled for 
the incident neutrino energy 

M 2 

E r u es = —2- = 4.2 x 10 21 eV (m^eV])" 1 , (17) 

1 1 It i 



7 



10 



26 



> 

CD 



10' 



CD 
CO 



LU 1( / 
LU 



25 



10 



AGASA 



Uniform sources 



io 19 io 20 
Energy [eV] 



Figure 4: Observed energy spectrum of cosmic rays extending beyond the GZK cutoff ffljj . 
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Figure 5: Energy spectrum of cosmic rays predicted from Z-bursts of extragalactic cosmic 



neutrinos and relic neutrinos with mass m v = 0.26 eV 221. 
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where Mz is the Z-boson mass. For neutrino masses 0(eV) these energies are remakably 
close to the highest energy cosmic rays observed. 



Z-bursts are indeed a possible explanation |2C], glj of the observed excess at energies 
above the Greisen-Zatsepin-Kuzmin (GZK) cutoff Eqzk = 4 x 10 19 GeV (fig. |). In a 
detailed study of extragalactic neutrinos a quantitive description of the data (fig. ^) is 



obtained for neutrino masses [22 



m v = 0.261&14 eV . (18) 

In fig. H the first bump at E = 4x 10 19 GeV represents protons produced at higher energies 
and accumulated just above the GZK cutoff. The second bum at E = 3 x 10 21 eV reflects 
the Z-burst. 

A confirmation of the Z-burst hypothesis would prove the existence of relic neutrinos 
and at the same time provide an absolute neutrino mass measurement! 



3 Weyl-, Dirac- and Majorana-Neutrinos 
3.1 C, P and CP 

The standard model is a chiral gauge theory, i.e. left- and right-handed fields have differ- 
ent electroweak interactions. Left-handed quarks and leptons are doublets with respect 
to the gauge group SU(2) of the weak interactions, whereas right-handed quarks and 
leptons are singlets, 





Ql = , , u R , d R , l L = , e R . (19) 



For massless particles helicity and chirality are identical. The chirality transformation 
corresponds to a multiplication of the Dirac spinor by 75, 

75^l(*,x) = -ip L (t,x) , 75?MM) = +ipR(t,x) . (20) 

Since left- and right-handed fields, which have different gauge interactions, also have 
opposite chiralities, the standard model is called a chiral gauge theory. 

For the standard model, the transformation properties of the fields under C, P and 
CP are of fundamental importance. Charge conjugation and parity are defined as 

C : %j}{t,x) ^^J c (t,x) = C^ T (t,x), (21) 
P : ip(t,x) ^ip p (t,x) = Pip(t,-x). (22) 
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Here C = 27270, ^ = V^7o and P = 70. For simplicity, we have set possible phase factors 
in the definition of C and P equal to I. One easily verifies that C and P transformations 
change chirality, 

7 5 ^f (t, x) = 757o^lO, -x) = +^f(t, x) , (23) 
75 ^(t, f) = 75^72^2 (*, x) = +^{t, x) . (24) 

Hence, parity and charge conjugated left-handed fermions are right-handed and vice 
versa. 

In general, a Dirac fermion can be decomposed into left- and right-handed parts, 

^ = ^ir^ + = i>r+i>l. (25) 

Parity and charge conjugation then relate these two components. However, in the stan- 
dard model left- and right-handed fields have different weak interactions. For the up- 
quark, for instance, one has, 

u = u R + u L , u L = ql , (26) 

i.e. ul, the left-handed partner of ur, changes under SU(2) gauge transformations. Hence, 
one cannot define P and C transformations which commute with the electroweak gauge 
symmetries. For neutrinos the situation is even worse, since right-handed neutrinos don't 
even exist in the minimal standard model! Only the CP transformation is well defined, 

CP : ip(t, x) -> ip CP (t, x) = CP^ T (t, -x) , (27) 

since it does not change chirality, 

TM# P = T*1% P = H>S P - (28) 

Hence, in principle, CP could be a symmetry of the standard model. 

Neutrinos are the quanta of the left-handed field vl(x). Neutrino states are created 
by the field operator which has the mode expansion, 

vl(x) = J dp(b(p)u L (p)e^ x + d\p)v L {p)e ip - x ) , (29) 

where dp = d 3 p/[(2ir) 3 2E}. The operators b\p) and ${p) create from the vacuum neu- 
trino and antineutrino states, respectively, 

\w,p) = tf(p)\0) , \v;p) = S(p)\0) ■ (30) 
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The operators for energy, momentum, helicity and lepton number are given by (cf. 



P, = -j d 3 x:^d^:= j dpp„ (tf(p)b(p) + S(jp)d(jp)) , (31) 
E • P = ^fd 3 x: :=- J dpE (b\p)b{p) - S{p)d{p)) , (32) 

L = I d 3 x : : = / dp (b\p)b{p) - d\p)d(p)) ; (33) 



here a 1 = ^ijkT '7 ■ Using the usual anticommutation relations for creation and annihi- 
lation operators one reads off from eqs. (|3~T|)-(|3"3"|), 

E • P 

p»W;p) =Pim\v;p) , —g-W',p) = -\f,p) = -i'\v,p), ( 34 ) 

E • P 

Pp\v;p) =p,j,\v;p) , — = +I^;p) = -£|^;p) , ( 35 ) 

i.e. the helicity (lepton number) of antineutrinos is positive (negative). Note, that helicity 
and lepton number always have opposite sign. 



3.2 Mass generation in the standard model 

Since the standard model is a chiral gauge theory quark and lepton masses can only be 
generated via spontaneous symmetry breaking. Starting from the Yukawa interactions 

Cy'^ = h ui jq Li URjHi + h di jq Li dRjH 2 + h ei jl Li e R jH2 + h.c. , (36) 

the vacuum expectation values of the Higgs fields, (H) 1 = v\ and (H) 2 = t>2, lead to the 
mass terms, 

£-M = m uijU L iU Rj + m dij d Li d Rj + m eij e Li e Rj + h.c. , (37) 
with the mass matrices for up-quarks, down-quarks and charged leptons, 

m u = h u v\ , m d = h d v 2 , m e = h e v 2 . (38) 
The mass matrices are diagonalized by bi-unitary transformations, 

V^m u V {u) = m dia9 , V {d ^m d V {d) = m d d iag , V {e)] m e V {e) = mf a9 , (39) 

with V^V^ = 1, etc. , which also define the transition from weak eigenstates ipLa, ipRa 
to mass eigenstates ipu, ip Ri , 

ULa = Vod U Li , = V^fd Li , . . . , e Ra = V^fe Ri . (40) 
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Since in general the transformation matrices are different for up- and down-quarks 
one obtains a mixing between mass eigenstates in the charged current (CC) weak inter- 
actions, 

4t = -JLY,UL a i»d La w; + ... 

= -^yisfVijdu w+ +..., (41) 

where Vy = V^Vjf is the familiar CKM mixing matrix. In general the Yukawa cou- 
plings are complex; hence, also the CKM matrix is complex, which leads to CP violation 
in weak interactions. 

Without right-handed neutrinos, weak and mass eigenstates can always be chosen 
to coincide for leptons. As a consequence, there is no mixing in the leptonic charged 
current, 

C% = ~\ E W""" W~ + ..., (42) 

* a 

and electron-, muon- and tau-number are separately conserved. On the contrary, in the 
quark sector only the total baryon number is conserved. 



3.3 Neutrino masses and mixings 

The simplest, and theoretically favoured, way to introduce neutrino masses makes use 
of right-handed neutrinos vr. This allows additional Yukawa couplings and a Majorana 
mass term, 

Cp = hvijluVRjHt + -Mij^VRj + h.c. (43) 

The Yukawa interaction defines the quantum numbers of the right-handed neutrino: it 
carries lepton number, which is a global charge, but no colour, weak isospin or hyper- 
charge, which are gauge quantum numbers. Hence, a Majorana mass term is allowed for 
the right-handed neutrinos, consistent with the gauge symmetries of the theory. It is very 
important that these masses are not generated by the Higgs mechanism and can there- 
fore be much larger than ordinay quark and lepton masses. This leads to light neutrino 
masses via the seesaw mechanism |[24|| . 

The Yukawa interaction which couples left-handed and right-handed neutrinos yields 
after spontaneous symmetry breaking the Dirac neutrino mass matrix mo = h u V\, so 
that the complete mass terms are given by 

Cm = m Dij VZiV Rj + -Mij^iyRj + h.c. (44) 
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In order to obtain the mass eigenstates one has to perform a unitary transformation. 
Using ITZmDVR = Vj l rnj D i/f j , the mass terms can be written in matrix form 



The unitary matrix which diagonalizes this mass matrix is easily constructed as power 
series in £ = mo/M. Up to terms 0(£ 3 ) one obtains 







["%) 





(46) 



In terms of the new left- and right-handed fields, L and R, the mass matrix is diagonal, 
£S? = i f I ) ( m : ° )( L :) + h,. (47) 



2 V J \ M I \ R 

Here the mass matrix m„ is given by 

m, = • (48) 

This is the famous seesaw mass relation p5]| . With M > m/), one obviously has 
m£). As an example, consider the case of just one generation. Choosing for mo the largest 
know fermion mass, mp ~ m t ~ 100 GeV, and for M the unification scale of unified 
theories, M ~ 10 15 GeV, one finds m u ~ 10~ 2 eV, which is precisely in the range of 
present experimental indications for neutrino masses. 

m u and M are the mass matrices of the light neutrinos and their heavy partners, 
respectively. The corresponding mass eigenstates are Majorana fermions, 

v = L + L c = v c , N = R c + R = N c . (49) 

L, R c and L c , R are the corresponding left- and right-handed components, 

L = ^—^v , R c = ^—t^N , etc. (50) 

_ _ 

In terms of the Majorana fields the mass terms read, 

£ { m = 2 Vmu ~Y^ u + 2 Nmv ~lP^ N + h ' c ' 
1 1— 

= -77(Re{m ; ,} + Tm{m i ,}75)z/ + -A r (Re{M} + Tm{M}75)A^. (51) 

_ _ 

Note, that in general the mass matrices have real and imaginary parts. This is a conse- 
quence of complex Yukawa couplings and a possible source of CP violation. 



13 



The mode expansion of the Majorana field operator reads, 



V (X) = V (X) 



dpJ2 (ai(p)ui(p)e- tp - x + al(p)vi(p)e ip - x 



(52) 



i=i 



where Ui(p) and Vi(p) = u^(p) are now solutions of the massive Dirac equation. In the 
massless case the connection with the left-handed neutrino field fl29|) is very simple, 



ai (p) = b{p) , a 2 (p) = d{p) , wi(p) = u L {p) , u 2 (p) = v c L {p) 



(53) 



i.e. the two spin states of the Majorana neutrino are identical with the neutrino and 
antineutrino states. In the massless case these states carry different lepton number, which 
is conserved. For massive Majorana neutrinos lepton number is violated by the Majorana 
mass term which couples the two polarization states. 

So far we have only discussed the simplest version of the seesaw mechanism. Ad- 



ditional Higgs fields also lead to a direct mass term for the left-handed neutrinos |26 



In principle, neutrino masses could also be generated by radiative corrections. This is a 
possible alternative to the seesaw mechanism. Particularly attractive are supersymmetric 



models with broken R-parity [p7[| , which lead to signatures testable at colliders. 

Neutrino masses imply mixing in the leptonic charged current. If the neutrino mass 
matrix is diagonalized by the unitary matrix U^ v \ 



( mi \ 

U^m v U (u) * = ~ m 2 

\ m 3 / 

the mixing matrix U appears in the charged current, 

.(0 _ 9 



(54) 



EW 



ij 



(55) 



where 



U - = U {e)] U {v) 

U «J u ia u aj 



(56) 



e-, fi- and r-number are now no longer separately conserved. The mixing matrix is 
frequently written in the following form, 



U 



I u el u e2 u e3 \ 

U,j,i U^2 U^s 

V u Tl u r2 u t3 J 



(57) 
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and can be expressed in the standard way as product of three rotation in the 12-, 13- 
and 23-planes, 



U 



/ 1 W 

C 23 S 2 3 
\ -S 23 C 23 / 



Cl3 





1 









Cl3 / 



\ / C12 S12 \ 
-Sl2 C12 

1 / 



y o o 



(58) 



As we shall see in the next section, present data are consistent with a small value of 
Sx3- The 3x3 mixing matrix then becomes a product of two 2x2 matrices describing 
mixing between the first and second, and the second and third generations, respectively. 
This very simple mixing pattern appears to be sufficient to account for the solar and the 
atmospheric neutrino deficits. 

4 Neutrino Oscillations 



In recent years there has been a wealth of experimental data in neutrino physics, and we 
can look forward to important new results also in the coming years. The present situation 
is summarized in fig. ^| which is taken from the review of particle physics. The latest most 
important result is the first data from SNO which we shall discuss in section 4.3. Detailed 



discussions of neutrino oscillations and references can be found in the reviews [29]- [32 . 



4.1 Oscillations in vacuum 



Neutrino oscillations |33], |34| are a very intriguing quantum mechanical effect similar 
to the well known oscillations between K°- and i^°-mesons. They occur because of the 
mixing in the charged weak current discussed in the previous section. The neutral and 
charged current weak interactions of neutrinos are described by the lagrangian 



C 



EW 



9 



9 - 



, , VLal^La Z ^ + -^=.e La ^U ai V Li W + h.C. 

2cosB w V2 



(59) 



where the fields eia, ot = 1 ... 3, represent the mass eigenstates of electron, muon and 
tau, and the fields z/ £i , i = 1 . . .n > 3, correspond to neutrino mass eigenstates. Hence, 
the charged lepton e a couples to the neutrino flavour eigenstate u a , which is a linear 
superposition of mass eigenstates, 



(60) 



Here v a and are spinors in flavour space, and U is a n x n unitary matrix. Three 
linear combinations of mass eigenstates have weak interactions, and are therefore called 
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Figure 6: The most important exclusion limits as well as preferred parameter regions 
from neutrino oscillation experiments assuming two-flavour oscillatons fftBj . 
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Figure 7: Neutrino production, propagation and absorption. 

active, whereas n — 3 linear combinations are sterile, i.e. they don't feel the weak force. 
In the case n — 4, for instance, the sterile neutrino is given by 



(61) 



In the following we will restrict ourselves to the case of three active neutrinos and shall 
not discuss the positive signature for oscillations from the LSND experiment J}£| at Los 
Alamos, which would require the existence of a sterile neutrino. 

Neutrinos are relativistic particles whose propagation is described by the Dirac equa- 
tion. Consider now a neutrino beam propagating in z-direction from a source at z = 
to a detector at z = L. The probability for neutrino oscillation is usually derived by 
considering first the time evolution of a spatially homogeneous state. Some hand waving 
arguments are then needed to obtain the physically relevant probability for oscillations 
in space, which, with lack of fortune, can also give the wrong result. As pointed out by 
Stodolsky, all this confusion can be avoided by realizing that the system corresponds to 
a stationary state |36], |37j which is completely characterized by the energy spectrum of 
the beam. For simplicity, we first consider the case of fixed energy E. The wave packet 



description of neutrino oscillations is discussed in ref. p8 



Neglecting the mixing between neutrinos and antineutrinos, which is suppressed by 
(rrii/E) 2 , the hamiltonian of the system reads in the mass eigenstate basis, 







/ \/p 2 + m\ 

+ 

V 



ml 








\ 











/ 



V 



2 

a mi 

o 







* ml 








a/j3 2 + m| J 



(62) 



P+i/ 

here p is the momentum operator conjugate to the z-coordinate, and we have assumed 
that the relevant eigenvalues are much larger than the neutrino masses raj. To describe 
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the neutrino oscillation u a — > z/^ we have to find an energy eigenstate, 

H u { E a \z) = EuP(z), (63) 
which satisfies the boundary condition 

4 a) (0) = ^ a . (64) 

The solution reads 



where the momentum eigenstates are given by 

= JPiZ 



v W(z)= e w Vil (66) 



with pi ~ E — mf/(2E). For the dependence of the wave function u% (z) on the z- 
coordinate one then obtains, 



E U ^ PiZ ^ ■ (67) 



Using the unitarity of the mixing matrix, WU = 1, this yields for the flavour transition 
probability between the source at z = and the detector at z = L, 

= ^ + ^W^(e- iAm ^-l), (68) 



where 



Aniy = — . (69) 

Terms with i = j do not contribute in the sum (|68"D , and a very useful form of transition 
probability is 



+ 2 ^ Im {U^UpiU^) sin 2 (Am^j . ( 70 ) 

Note the different oscillation lengths of the real and imaginary parts. In the realistic case 
of a neutrino beam with some energy spectrum p(E), one has to integrate over energy, 

P(v a -+ Vf} ) = J dEp{E)P{v a -> up) . (71) 

All patterns of neutrino oscillations including CP violation are described by the 
master formula (|70"|). Let us discuss some of its properties: 
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The dependence on L shows the expected oscillatory behaviour, which is an inter- 
ference effect and dissapears for Am,?- = 0. 

Due to the unitarity of the mixing matrix U the total flux of neutrinos is conserved, 

(3 i 

In the simplest case of only 2 flavours (a, (3 = e, fi) only one term contributes in 
the sum of (|70|) . The unitarity of the mixing matrix U then implies (i = 2,j — 1), 
Ua2Up2U a iUg 1 = — | U a 2\ 2 \U/S2 \ 2 ■ This yields for the transition probability 



P (y a ^ Up ) = 4\U a2 \ A \U 02 \ 2 sm 2 (Am 2 — J , (73) 

where Am 2 = m 2 , — m\. The unitary 2x2 mixing matrix can be parametrized in 
the standard way in terms of 1 rotation angle and 3 phases, 

_ / e^ 1 cos O e^* 2- * 3 ) sin O \ 
' I -e 4W,1+<fe) sin6o e^ 2 cos6 J 

With 4 1 ZTq-2 1 2 1 1 2 = sin 2 2Bo one then reads off from eq. (|73|) the standard formula 
for the transition probability in the case of 2 flavours, 

P(u a -> up) = sin 2 29 sin 2 — - , (75) 

Llr. 



■'vac 



with the vacuum oscillation length 

AE 



Am? Ap 

E[GeV] 
Am 2 [eV 2 



1-27 km hffljT • ( 76 ) 



For fast oscillations, i.e. large oscillation phase, averaging over the energy resolution 
of the detector and the distance L according to the uncertainty of the production 
point, yields the average transition probability, 

1 



P{v a ^v p ) = -sin 2 26 . (77) 

In disappearance experiments one measures the survival probability P(u a — > u a ) 
which is directly related to the transition probability, 

P(u a -+v a ) = 1 - P(v a -> Up) 

= 1 - A\U a2 \ 2 (l - \U a2 \ 2 ) sin 2 ^Am^j . (78) 
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L[km] 


E[GeV] 


Am 2 [cV 2 ] 


accelerator (short baseline) 


0.1 


1 


10 


reactor 


0.1 


lO- 3 


io- 2 


accelerator (long baseline) 


10 3 


10 


10~ 2 


atmospheric 


10 4 


1 


10~ 4 


solar 


10 8 


io- 3 


io- 11 



Table 1: The approximate reach in Am 2 of different oscillation experiments. 

Particularly interesting is the case of three hierarchical neutrinos, m\ -C m 2 <C m 2 . 
Suppose that Am 2 1 L/(4S) = 0(1), and correspondingly Am 2 1 L/(4£') <C 1. From 
eq. (|70| ) one then obtains 

P(u a ^u p ) ~ -4Re(f/: 3 ^3^2^2 + ^3^3^i^i) sin2 ( Am 3iT|) 

= 4|^ 3 ri^ 3 rsin 2 ^Am 2 1 ^ . (79) 

Note, that this result is completely analogous to the two-flavour case; only the large 
mass difference matters. 

The sensitivity of an oscillation experiment with respect to the neutrino mass 
difference Am 2 is determined by the neutrino energy E and the oscillation length 
L. Oscillations become visible for mass differences above Am 2 L/(4E) ~ 1. In the 
appropriate units of energy and length this condition reads 

Am 2 [eV 2 ] . |M . m 

Relevant examples are given in table |T|. 

In the past the sensitivity of neutrino oscillation experiments has been mostly 
displayed in the (Am 2 , sin 2 26) plane. Obviously, this parametrization covers only 
half of the parameter space. More appropriate are the variables (Am 2 , tan 2 ©), 



which is particularly important for oscillations in matter p9 . 



The measurement of magnitude and sign of CP violation in neutrino oscillations 
would be of fundamental importance. One CP violating observable is the asymme- 
try 

A eiM = P{u e -> u„) - P(y e -> z/ M ) . (81) 
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Using CPT invariance and the unitarity of the mixing matrix one easily verifies 
that in the case of three neutrino flavours the three possible asymmetries are all 
equal, 

A efl = A MT = A re . (82) 
From the master formula (|70|) one reads off, 

A Q/3 = 4 J2 Im (KPfiUvjUpi) sin2 ( Am ^) • ( 83 ) 

i>j ^ ' 

The product of mixing matrices is completely antisymmetric in the flavour indices, 

Im (JTJJpiUajU},) = eaf&jJi , (84) 

where ? a b = Ylc=i e abc and Ji is the leptonic Jarlskog parameter. For quarks one 
finds for the corresponding quantity J q ~ 10~ 5 . Due to the large neutrino mixings 
Ji turns out to be much larger, which gives hope to observe CP violating effects in 
future superbeam experiments and at a neutrino factory GO . 



4.2 Oscillations in matter 

In vacuum, neutrino oscillation probabilities are bounded by the mixing angle, P < 
sin 2 20 o . Hence, for small mixing angles transition probabilites are small. In matter, a 
resonance enhancement of neutrino oscillations can take place and transition probabilities 
can be maximal even for small vacuum mixing angles — this is the Mikheyev-Smirnov- 



Wolfenstein effect H 



The matter enhancement of oscillations is a coherent effect, due to elastic forward 
scattering of neutrinos with negligible momentum transfer. The effect can be described 
by the propagation of the neutrinos in an approximately constant potential generated 
by the exchange of H^-bosons (fig. §), 



->cc 



-2V2G F eE^v Le 17Z^e L . (85) 



The corresponding Z-exchange does not distinguish between i/ e , and u T , and therefore 
it has no influence on neutrino oscillations. 



After a Fierz-transformation one obtains from eq. fl8~5|) , 

L C C = -2\f2G F VTel il VL<^Ll l L e L ■ ( 86 ) 



In order to obtain the effective potential for the neutrino propagation in matter, e.g. the 
interior of the sun, one has to evaluate the expectation value of the electron current in 
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Figure 8: Neutrino interaction with matter via charged and neutral currents. 

the corresponding state, i.e. the quantity (6^(1— j 5 )e). Since no polarization vector and 
no spatial direction are singled out, one has 

(e 7 ^ 75 e) = , (eYe) = , (87) 

(e 7 °e) = ( e t e ) = 2n e , (88) 

where n e is the electron number density and we have summed over electron spins. 

The final lagrangian describing the propagation of neutrinos now reads (cf. ([51]), 

m v = ml), 

1 Gp 

For simplicity we shall restrict ourselves in the following to the case of two flavours, i.e. 
a, (3 = e,/z. We also consider the idealized case of constant electron density. From ( |89| ) 
one obtains the hamiltonian 

H = H + V , (90) 

with 

U 0^1^ rr ( V2G F n e 0\ 

H = 7 7-<9 + i rn , V = . 91 

% \ J 

Note, that the matter induced potential has the same Lorentz structure as a static 
Coulomb potential. Hence, matter effects have opposite sign for neutrinos and anti- 
neutrinos. 

The probability for flavour oscillations can now be calculated completely analogous 
to the case of vacuum oscillations. We again have to find a stationary state, 

Hvf{z)=Evf{z), (92) 
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which satisfies the boundary condition, 

4 a) (o) = ^. 

The free hamiltonian is known in the mass eigenstate basis (cf. where we again 

neglect the mixing left- and right-handed states. In the flavour basis one then has 

H = U(<d )H U\e ) + V . (93) 

With 



U(Q )=( C0S9 ° Sin6 ° 
y — sin 6 cos 6 



(94) 



/ p + - ^ cos 29 + V2G F n e ^ sin 29 \ 

Am 2 • ,n - , ml+ml Am 2 mo oa " ^ 



the hamiltonian reads explicitely 

-^cos 

^sin26 p+%^ + ^cos20 o 

Here Am 2 = m 2 — m|, and we have assumed that the relevant eigenvalues of p are 
E + O (mf/E). It is now straightforward to diagonalize H, requiring 

W(Q)HU(Q) =( E Q ° J . (96) 



From (BH) one obtains the rotation angle 



^sin26o , N 

tan29 = - r ^ — ^ , 97 

^cos2e - V2G F n e 



and the equation for the momentum eigenvalues 
Pi , 1£ - 



+ + >/2Cfrn. - # 



P2 + ^ + V2G F n e + | 



4£ 

where 



)-l 








v o 





(98) 



/\ ?T?,^ 

A = — — cos 20 o - V2G F n e . (99) 
2£< 

This yields the momentum eigenvalues 

Pl +p 2 = 2E - "^"y"' ~ V2G F n e , (100) 

(/Am 2 \ 2 /Am 2 \ 2 \ 1//2 

( -^cos2Q -V2G F n e ) +(^ sin20 o) J -(101) 
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Figure 9: Change of neutrino momentum and flavour in the sun. 



The probability for flavour oscillations is again given by eq. (|75|) , with Go replaced 
by 6 and L vac replaced by L mat = 2/Ap, 

~ sin 2 26 sin 2 — ^— . (102) 

Lmat 

The expression ( |9~7D for the rotation angle has the typical resonance form, and maximal 
mixing, i.e. = 45°, is reached for 

^2G F n e = — — cos26 . (103) 
2E 

This is the M SW resonance condition. It requires in particular that cos 2Go is posi- 
tive. A negative sign would allow resonant oscillations for antineutrinos. 

At the center of the sun the electron density is n e ~ 10 26 /cm 3 . With E ~ 1 MeV, 
the MSW resonance condition can be fullfilled on the way out to the surface of the sun 
for mass differences Am 2 ~ 10 -5 eV 2 . In the 'adiabatic approximation', where the oscil- 
lations are fast compared to the variation of the electron density, eqs. ( |100|) , (|1Q1|) and 



( |103| ) can be used at each distance z from the center. The physical picture is particularly 
simple in the case of small mixing angle Go, where the off-diagonal terms in the hamil- 
tonian fl95|) are small. The momentum of the produced electron neutrino increases, and 
near zr the transition v e — > takes place after which the momentum stays constant 
(fig. |]). A more accurate calculation of the oscillation probability can be carried out 
similar to the treatment of the Landau-Zener effect in atomic physics []30| , [32|. 
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Figure 10: Solar neutrino spectrum and energy thresholds of solar neutrino experiments 



4.3 Comparison with experiment 

We are now ready to interpret the deficit in the solar and atmospheric neutrino fluxes in 
terms of neutrino oscillations. In the following only a brief account of the main results 
will be given. More detailed discussions of the fusion processes in the sun, solar neutrino 
experiments, the atmospheric neutrino anomaly, reactor experiments and global fits can 
be found in [H]-[|32|. 



4.3.1 Solar neutrinos 

The solar energy is generated by thermonuclear reactions, in particular the pp cycle 
and the CNO cycle, which also produce electron neutrinos. The largest neutrino flux, 
$i, e = 6.0 x 10 10 /[cm 2 s] is due to the reaction, 



p + p^d + e + + v e 



(104) 



with a maximum neutrino energy of 0.42 MeV. Also important are the processes e + 
7 Be ^u e + 7 Li with $„ e = 4.9 x 10 9 /[cm 2 s], E v = 0.86 MeV and 8 B -> 8 Be + e+ + u e , 
with $j, e = 5.0 x 10 6 /[cm 2 s], E v < 0.86 MeV. The complete energy spectrum of solar 
neutrinos is summarized in fig. HI] . 



The deficit in the solar neutrino flux was discovered in radiochemical experiments. 
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Total Rates: Standard Model vs. Experiment 
Bahcall-Pinsonneault 2000 




Figure 11: Theoretical flux predictions compared with experimental results; for the CI and 
Ga experiments the units are SNU, for Kamiokande and Superkamiokande the ratio of 



measured and predicted fluxes are shown fcffi] . 



The first one was the Homestake experiment of Davis et al. based on the reaction 

z/ e + 37 Cl^ 37 Ar + e" , (105) 

with an energy threshold of 0.814 MeV. Only 1/3 of the predicted flux was observed. 
The GALLEX, SAGE and GNO experiments detect solar neutrinos via the reaction 

v e + 71 Ga -> 71 Ge + e~ . (106) 

Since the energy threshold is only 0.234 MeV, these experiments could see for the first 
time neutrinos from the dominant process ( |104|) . About 60% of the predicted flux is 
observed in this reaction (fig. |TT|). 

The Kamiokande and Super-Kamiokande experiments are water Cherenkov detec- 
tors. Here the reaction, 

v% + e~ -»• u x + e~ , (107) 

is used to detect solar neutrinos. To suppress background only high-energy neutrinos 
from the 8 B decay, with E > 7 GeV, are detected. As a consequence the recoil electrons 
are peaked in the direction of the incoming neutrino. The angular distribution of the 
detected electrons shows indeed a peak opposite to the direction to the sun. This clearly 
demonstrates the solar origin of detected neutrinos! 
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Figure 12: Elastic neutrino electron scattering via charged and neutral currrents. 



An important step towards the final resolution of the solar neutrino problem was 
made by the recent results of the Sudbury Neutrino Observatory (SNO). Like in Super- 
Kamiokande, high-energy S B neutrinos are observed, this time in a 1000 t detector of 
heavy water D 2 0. This allows the observation of the charged current (CC) reaction, 



v e + d^p + p + e~ 
as well as the elastic scattering (ES), 

v a + e~ -> Va + e~ 



(108) 



(109) 



which involves charged and neutral currents (fig. 12). From the CC reaction ( |108|) one 
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Figure 13: Two-flavour (v e to u a ) fit of solar neutrino data, with confidence levels 90%, 
95%, 99%, 99.7% 0/. 
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v s fits, leaving the 8 B flux free 



can, for the first time (!), directly determine the v e flux [}43 



$ cc (z/ e ) = 1.75 ± 0.07 (stat^o.nisy st) ± 0M(theo) x lO^cm" 



\+0.12 



(110) 



Here statistical, systematic and theoretical errors are listed separately. Under the as- 
sumption that the entire flux consists of electron neutrinos, one obtains from the ES 
process [ |4"3f . 

§ ES (v e ) = 2.39 ± OM(stat)t it{syst) x lf/cm-V 1 . (Ill) 

Comparing Q cc and Q ES one concludes that with a significance of 1.6 a there is evidence 
for a non-f e component in the solar neutrino flux. Combining this with the more precise 
measurement of $ by Superkamiokande the significance increases to 3.3 a. 

Alternatively, one may assume that electron-neutrinos have partially been converted 
into muon- and tau-neutrinos. In this case, the flux extracted from the elastic scattering 
( |109| ) should coincide with the flux of 8 B neutrinos. One finds [43j|, 



5.44 ± 0.99 x 10 6 cm- 2 s _1 



(112) 



in remarkable agreement with the prediction of the standard solar model! 

Global fits of all solar neutrino data have been performed under different assumptions 



44, 45 1 . Fig. 13 shows a two-flavour fit for the oscillation of v e into an active neutrino v a 



using the solar flux PB2000 (fig. [To]) . Four regions in parameter can describe the data: 
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LMA (Large Mixing Angle MSW solution), SMA (Small Mixing Angle MSW solution), 
LOW (MSW solution with low Am 2 ) and VAC (Vacuum oscillation solution). The LMA 
solution gives the best fit. 

In fig. |14] two-flavour fits for oscillations into active and sterile neutrinos are com- 
pared. Clearly, the oscillation v e —>■ v s is essentially excluded. For the oscillation into a 
superposition of active and sterile neutrino the best fit is obtained for vanishing sterile 
component. 



4.3.2 Atmospheric neutrinos 

Cosmic rays produce in the earth's atmosphere pions and kaons which decay into charged 
leptons and neutrinos, 

M ± -> e ± + v e (v e ) + Vptyn) . (113) 

These neutrinos can be detected by the usual neutrino-nucleon charged current reactions. 
Naive counting suggests that the corresponding flux of 'atmospheric neutrinos' contains 
twice as many muon-neutrinos than electron-neutrinos. More accurately, one compares 
the measured ratio of v^s and z/ e 's, N^/Ne, with the theoretical expectation of this 
quantity. For the double ratio, 

p _ (N^/N e ) data (-\-\A\ 



the Super-Kamiokande collaboration finds 

R = 0.638 ± 0.017(stat) ± 0.050(syst) , (115) 

where the sample has been restricted to sub-GeV (E <~ 1 GeV) charged leptons. For 
the sample of multi-GeV (E >~ 1 GeV) leptons a similar number is found. 

It appears natural to interpret the observed deficit again in terms of neutrino 
oscillations. This hypothesis is further strengthened by the analysis of the zenith angle 
dependence shown in fig. [15]. No deficit is seen for electron-neutrinos. On the contrary, 
for muon-neutrinos the deficit is the larger, the larger the zenith angle, i.e. the larger the 
distance between detector and production point in the earth's atmosphere. For upward 
going neutrinos the oscillation length is about 13.000 km. 

For neutrinos with energies O(GeV) one has a sensitivity for mass splittings down 
to Am 2 ~ 10~ 4 eV 2 . A two-flavour analysis — > v a yields for the allowed parameter 
range at 90% CL, 

sin 2 26 > 0.88 , 1.6 x 10~ 3 < Am 2 < 4 x 10~ 3 . (116) 
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Figure 15: Zenith angle distribution of Super- Kamiokande 1289 day sub-GeV and multi- 
GeV samples. Solid and dashed lines correspond to MC with no oscillation and MC with 
best fit oscillation parameters, respectively 
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For oscillations into sterile neutrinos, — > u s , earth matter effects become important. 
A detailed analysis shows that this case is disfavoured. 

Important information on neutrino mixing is also obtained from reactor experiments. 



Failure to observe v e — ► u x oscillations gives the bound [[17]], for Ami ~ 3 x 10 3 eV 2 , 



sin 2 26 < 0.12 . (117) 

For three active neutrinos, this suggests a small mixing between the first and the third 
generation. 



4.3.3 Future prospects 

In the coming years we can look forward to important results from several experiments. 
Mini-BooNE at Fermilab will verify or falsify the LSND signal for oscillations. The long- 
baseline experiments K2K and MINOS will be able to verify — v a oscillations and 
to determine the mass difference Am 2 3 more precisely. Correspondingly, OPERA and 
ICARUS at Gran Sasso, using the muon-neutrino beam from CERN, will look for r 
appearance, verifying unequivocally — v T oscillations. The first oscillation dip could 
be studied by MONOLITH. The 1 Be line of solar neutrinos will be studied by Borexino. 

KamLAND is the first terrestrial experiment sensitive to the solar neutrino signal 
(fig. [Tj^). Using reactor neutrinos over a baseline of about 175 km it could verify the 
LMA solution of the solar neutrino problem. This would be of crucial importance for the 
prospects of observing CP violation in neutrino oscillations. This may then be possible 
in neutrino superbeam experiments (JHF in Japan, SPL at CERN) or, eventually, at a 
neutrino factory. 

The result of neutrino oscillation experiments will be the determination of neutrino 
mass differences and the leptonic mixing matrix. At present, a three-flavour analysis 
assuming the LMA solution leads to the following result p9| , 



U 



( 0.74- 0.90 0.45- 0.65 < 0.16 \ 

0.22- 0.61 0.46- 0.77 0.57- 0.71 
\ 0.14 - 0.55 0.36- 0.68 0.71 - 0.82 / 



(118) 



The emerging structure is very remarkable and completely different from the CKM matrix 
of quark mixing. Whereas Vqkm is strongly hierarchical, all elements of the leptonic 
mixing matrix U are of the same order, except U e s, for which only an upper bound 
exists. This has important implications for the structure of fermion masses in unified 
theories. 



31 




Figure 16: Sensitivity of the KamLAND neutrino oscillation experiment. The curves 
represent 90% CL limits which can be reached in 3 years of running. The cases a) - c) 
correspond to different assumptions about the background 

5 Neutrino Masses in GUTs 



Neutrino masses and their relation to quark and charged lepton masses play an impor- 
tant role in grand unified theories. According to the seesaw mechanism the smallness of 
the light neutrino masses is explained by the largeness of the heavy Majorana neutrino 
masses, which extend up to the mass scale of unification. Hence, the present experimen- 
tal evidence for neutrino masses and mixings probes for the first time the physics of 
unification. Reviews and extensive references are given in ref. |5U||, a special approach is 



described in ref. |jT|. 



5.1 Elements of grand unified theories 

Let us first consider the basic ingredients of GUTs, starting from the symmetries of 
the standard model where neutrinos are massless. As a consequence, four 'charges' are 
classically conserved, three lepton numbers and baryon number, 

L e , L M , L T , B. (119) 

The Adler-Bell-Jackiw anomaly (fig. [H]), a quantum effect, reduces the four conserved 
charges to three [52|, which may be chosen as, 

L e — L^, L^ — L T , B — L, (120) 
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Figure 17: A B J- anomaly of the B — L current. 

where L = L e + L M + L T is the total lepton number. One may wonder whether these 
charges correspond to fundamental global symmetries, or whether they are just acciden- 
tal symmetries of the low-energy effective theory, similar to weak isospin which is an 
approximate symmetry of weak interactions at energies much below the W-boson mass. 

Neutrino masses and mixings break the first two of these symmetries, L e — and 
— L T , but presently we do not know whether B — L is a fundamental symmetry or 

not. B — L is conserved if neutrinos have only Dirac masses and the Major ana masses of 

the right-handed neutrinos vanish, 

m D = h v {H x ) , M = . (121) 

Note, that this requires tiny Yukawa couplings, e.g. h u2 2 ~ 1CT 14 for m V 2 ~5x 10~ 3 eV. 
This possibility can presently not be excluded although it might appear unnatural. Al- 
ternatively, B — L may be broken at some intermediate mass scale or at the GUT scale 
A ~ 10 16 GeV, where the gauge couplings unify in the supersymmetric standard model 



(fig. [Lq) . This is indeed suggested by the smallness of the observed neutrino masses in 
connection with the seesaw mechanism. With h v ?& = 0(1) one obtains for the mass of 
the heaviest Majorana neutrino from eq. 



ITT \2 

M 3 ~ ~ 10 15 GeV ~ Agut ■ (122) 



The unification of gauge couplings suggests that the standard model gauge group is 
part of a larger simple group, 

Gsm = U(l) x SU(2) x SU(3) C SU(5) C . . . . (123) 

The simplest GUT is based on the gauge group SU(5) p3| . Here quarks and leptons are 
grouped into the multiplets, 

10 = (q L , u% e c R ) , 5* = (d c R , l L ), l = u R . (124) 
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Figure 18: Unification of gauge couplings in the supersymmetric standard model. 

Note that, unlike gauge fields, quarks and leptons are not unified in a single irreducible 
representation. In particular, the right-handed neutrinos are gauge singlets and can there- 
fore have Majorana masses not generated by spontaneous symmetry breaking. In addition 
one has three Yukawa interactions, which couple the fermions to the Higgs fields -Hi (5) 
and F 2 (5*), 

£ = /w, 1040^1(5) + ^5**10^2(5*) + K^djH^b) + MyUh • (125) 

The mass matrices of up-quarks, down-quarks, charged leptons and the Dirac neutrino 
mass matrix are given by m u = h u V\, nid = h&Vi = m e and mo — h v v\, respectively, 
with v\ = (Hi) and v 2 = (H 2 ). The SU(5) mass relation = m e is successful for the 
third generation but requires substantial corrections for the second and first generations. 
The Majorana masses M are independent of the Higgs mechanism and can therefore be 
much larger than the electroweak scale v. 

Once right-handed neutrinos are introduced, the B — L charges of each generation 
add up to zero. Hence, there are no mixed B — L — gravitational anomalies (fig. O), and 
the U(1)b-l symmetry can be embedded together with SU(5) into a larger GUT group. 

In this way one arrives at the gauge group SO(10) [j54|j . All quarks and leptons of 
one generation are now contained in a single multiplet, 

16 = (q L , u% e c R , d c R , l L , v R ) . (126) 

Quark and lepton mass matrices are obtained from the couplings of the fermion multiplets 
16i to the Higgs multiplets #i(10), if 2 (10) and $(126), 

C = KijlGilGjHxilQ) + ^1646^2(10) + /i iV iil6 l 16 i $(126) . (127) 
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Figure 19: Mixed B — L — gravitational anomaly. 



Here we have assumed that the two Higgs doublets of the standard model are contained in 
the two 10-plets Hi and H 2 , respectively. This yields the quark mass matrices m u = h u V\, 
Tn>d = hdV2, with v\ = (Hi) and v 2 = (H 2 ), and the lepton mass matrices 



m,. 



VTLf, 



m d . 



(128) 



Contrary to SU(5) GUTs, the Dirac neutrino and the up-quark mass matrices are now 
related. Note, that all matrices are symmetric. The Majorana mass matrix M = Jin ($), 
which is also generated by spontaneous symmetry, is a priori independent of m u and m^. 

In the literature also GUT groups larger than SO(10) have been studied. Particularly 
interesting is the sequence of exceptional groups which terminates at rank 8, 



£ 4 = SU(5) CE 5 = SO (10) CE 6 CE 7 CE 8 . 



(129) 



The last two groups can also unify different generations and thereby restrict the Yukawa 
matrices. They arise naturally in higher dimensional supergravity theories and in string 
theories. 

What are the GUT predictions for neutrino masses and mixings? As emphasized at 
the end of section 4.3.3, the emerging structure of the leptonic mixing matrix appears 
to be remarkably simple, 

/ * * o \ 

* * * 

* * * 



U 



\ 



(130) 



/ 



where the V denotes matrix elements whose value is consistent with the range 0.5 . . . 0.8, 
whereas for the matrix element 'o' only an upper bound exits, \U e s\ < 0.16. 

There are several interesting 'Ansatze' to explain this pattern, such as 'bi-maximal' 
mixing (0i2 = 6 23 = 45°, 6 13 = 0), 
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or 'democratic mixing' where in the weak eigenstate basis all elements are equal to 1, 
which leads to the mixing matrix, 
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Further, one can also construct 'tri-bimaximal' mixing [|55]. All these patterns are inter- 
esting. In general, however, they lack an underlying symmetry. In the following we shall 
restrict our discussion to constraints on neutrino masses which arise for the simplest 
GUT groups, SU(5) and SO(10). 

An important consequence of neutrino mixing are flavour changing processes and 



electric dipole moments of charged leptons j56]. In particular supersymmetric theories 
predict effects large enough to be discovered in the near future, even before the start of 
LHC. 



5.2 Models with SU(5) 

An attractive framework to explain the observed mass hierarchies of quarks and 
charged leptons is the Froggatt-Nielsen mechanism ]57 based on a spontaneously bro- 



ken U(l)j? generation symmetry. The Yukawa couplings are assumed to arise from non- 
renormalizable interactions after a gauge singlet field $ acquires a vacuum expectation 
value, 

f ) • (133) 

Here g^ are couplings 0(1), and Qi are the U(l)j? charges of the various fermions, with 
Q$ = — 1- The interaction scale A is usually chosen to be very large, A > A GUT . 

The symmetry group SU(5) xU(l) F has been considered by a number of authors. 
Particularly interesting is the case with a 'lopsided' family structure where the chiral 
U(l)i? charges are different for the 5*-plets and the 10-plets of the same family f58|- 



folf . Note, that such lopsided charge assignments are not consistent with the embedding 
into a higher-dimensional gauge group, like SO(10)xU(l)^ or E s xU(l)p. An example of 
phenomenologically allowed lopsided charges Qi is given in table |2|. 
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Table 2: Lopsided U(l)p charges of SU(5) multiplets. 
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The charge assignements determine the structure of the Yukawa matrices, e.g., 



h e , h v ~ 



e 3 e 2 e 2 \ 



e e e 
ell/ 



(134) 



where the parameter e = ($) /A controls the flavour mixing, and coefficients 0(1) are 
unknown. The corresponding mass hierarchies for up-quarks, down-quarks and charged 
leptons are, 



m t : m c : m u ~ 
m b :m s :m d = 



1 : e 2 : e 4 



1 : e : e d 



(135) 
(136) 



The differences between the observed down-quark mass hierarchy and the charged lepton 



mass hierarchy can be accounted for by introducing additional Higgs fields [ 62] . From a fit 
to the running quark and lepton masses at the GUT scale the flavour mixing parameter 
is determined as e ~ 0.06. 

The light neutrino mass matrix is obtained from the seesaw formula, 
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Note, that the structure of this matrix is determined by the U(l)j? charges of the 5*-plets 
only. It is independent of the \J(1)f charges of the right-handed neutrinos. 



Since all elements of the 2-3 submatrix of ( 137 ) are 0(1), one naturally obtains a 
large — v T mixing angle 9 2 3 0, ^Pfl . At first sight one may expect G>i 2 = 0(e), which 
would correspond to the SMA solution of the MSW effect. However, one can also have 



a large mixing angle G12 if the determinant of the 2-3 submatrix of m v is 0(e) |63 



Choosing the coefficients 0(1) randomly, in the spirit of 'flavour anarchy' [p9|| , the SMA 
and the LMA solutions are about equally probable for e ~ 0.1 p3. The corresponding 
neutrino masses are consistent with m 2 ~ 5 x 10~ 3 eV and m 3 ~ 5 x 10~ 2 eV. We 
conclude that the neutrino mass matrix ( |137| ) naturally yields a large angle 623, with 
012 being either large or small. In order to have maximal mixings the coefficients 0(1) 
have to obey special relations. 



5.3 Models with SO(10) 



Neutrino masses in SO(10) GUTs have been extensively discussed in the literature |65 . 
Since all quarks and leptons are unified in a single multiplet these models often have 
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difficulties to reconcile the large neutrino mixings with the small quark mixings. In the 
following we shall illustrate this problem by means of an example [^BJ which only makes 
use of the seesaw relation, the SO(10) relation between the up-quark and Dirac neutrino 
mass matrices, 



m u = m D , 

and the empirically known properties of the up-quark mass matrix. 
With m d = m u the seesaw mass relation becomes 



m v ~ -m u —m u 



(138) 



(139) 



The large neutrino mixings now appear very puzzling, since the quark mass matrices are 
hierarchical and the quark mixings are small. It turns out, however, that because of the 
known properties of the up-quark mass matrix this puzzle can be resolved provided the 
heavy neutrino masses also obey a specific hierarchy. This then leads to predictions for a 
number of observables in neutrino physics including the cosmological baryon asymmetry. 

From the phenomenology of weak decays we know that the quark matrices have 
approximately the form ||67|| , 
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Here e 1 is the parameter which determines the flavour mixing, and p = \p\e ta , 
7] = \r)\e l P are complex parameters 0(1). We have chosen a 'hierarchical' basis, where 
off-diagonal matrix elements are small compared to the product of the corresponding 



eigenvalues, 



m. 



i.i i 



< 0(\rriimj\). In contrast to the usual assumption of hermitian mass 



matrices, SO(10) invariance dictates the matrices to be symmetric. All parameters may 
take different values for up— and down— quarks. Typical choices for e are e u ~ 0.07, 
€d ~ 0.2. The agreement with data can be improved by adding in the 1-3 element a term 
0(e A ) which, however, is not important for the following analysis. In case of a 1-3 element 
0(e 3 ), one can have a SU(3) generation symmetry leading to different results pSj. Data 
also fix one product of phases to be 'maximal', i.e. A = <p u — a u — (f>d + a<2 — n/2. 

We do not know the structure of the Majorana mass matrix M = ($). However, 
in models with family symmetries it should be similar to the quark mass matrices, i.e. 
the structure should be independent of the Higgs field. In this case, one expects 
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with M12 < M 22 
unitary matrix, 



M23 <C M33. The symmetric mass matrix M is diagonalized by a 
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M 3 / 



(142) 



Using the seesaw formula one can now evaluate the light neutrino mass matrix. Since 
the choice of the Majorana matrix M fixes a basis for the right-handed neutrinos the 
allowed phase redefinitions of the Dirac mass matrix mp are restricted. In eq. (|140|) the 
phases of all matrix elements have therefore been kept. 

The v il -v T mixing angle is known to be large. This leads us to require m Vi . = 0(1) 
for i, j = 2, 3. It is remarkable that this determines the hierarchy of the heavy Majorana 
mass matrix to be 



M 12 : M„ : M 
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'33 



e 4 :l 



With M 33 ~ M 3 , M 22 = <xe 4 M 3 , M 2 
masses and mixings to order 0(e 4 ), 



Ce 4 M 3 ~ M 22 and M 12 



(143) 



e 5 Ms, one obtains for 



Mi ~ M 3 

a 



Mo ~ cre 4 M3 



(144) 
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(9(e 4 ) and U- 
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0. Note, that a can always be 



chosen real. This yields for the light neutrino mass matrix 





( 





ee 2 ^ 





\ 


m v = - 




ee 2i<P 


-ae 2i<p + 2pe i * 








\ 





rje 1 ^ 




J 



± 



(145) 



The complex parameter ( does not enter because of the hierarchy. The matrix (|145|) 
has the same structure as the mass matrix ( |137| ) in the SU(5)xU(1)f model, except for 
additional texture zeroes. Since, as required, all elements of the 2-3 submatrix are 0(1), 
the mixing angle 023 is naturally large. A large mixing angle B12 can again occur in case 
of a small determinant of the 2-3 submatrix. Such a condition can be fullfilled without 
fine tuning for a, p,rj = 0(1). 

The mass matrix m v is again diagonalized by a unitary matrix, U^m v U^* = 
diag(mi, m 2 , m 3 ). A straightforward calculation yields (s^ = sinOjj, Cy = cosGy, £ = 

^/(1 + M 2 )), 
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with the mixing angles, 

tan20 23 ~ - ^ .„ , tan26 12 ~ 2^1 + |r/| 2 ^ . (147) 
1 — |r]| z o 

Note, that the 1-3 element of the mixing matrix is small, U{£ = 0(e). The masses of the 
light neutrinos are 

6 

mi~- tan 2 G>i 2 m 2 , m 2 ~ , cot 9 12 m 3 , m 3 ~ (1 + |??| 2 ) — . (148) 

(1 + \r]\ z ) 6 ' z M 3 

This corresponds to the weak hierarchy, 

mi : m 2 : m 3 = e : e : 1 , (149) 

with m 2 ~ m 2 ~ Am^ = m 2 — m 2 ~ e 2 . Since e ~ 0.1, this pattern is consistent with 
the LMA solution of the solar neutrino problem, but not with the LOW solution. 

The large z/ M -z/ r mixing is related to the very large mass hierarchy ( |144j ) of the heavy 
Majorana neutrinos. The large v^-v^ mixing follows from the particular values of parame- 
ters 0(1). Hence, one expects two large mixing angles, but single maximal or bi-maximal 
mixing would require fine tuning. On the other hand, one definite prediction is the oc- 
curence of exactly one small matrix element, u[ 3 = 0(e). Note, that the obtained pattern 
of neutrino mixings is independent of the off-diagonal elements of the mass matrix M. 
For instance, replacing the texture ( |141| ) by a diagonal matrix, M = diag(M l5 M 2 , M 3 ), 
leads to the same pattern of neutrino mixings. 

In order to calculate various observables in neutrino physics we need the leptonic 
mixing matrix 

U = U {e)] U {v) , (150) 

where is the charged lepton mixing matrix. In our framework we expect ~ V®, 
and also V = V^V^ ~ for the CKM matrix since e u < e<f. This yields for the 
leptonic mixing matrix 

U ~ V ] U {U) . (151) 

To leading order in the Cabibbo angle A ~ 0.2 we only need the off-diagonal elements 
= A = — V 2 \ ■ Since the matrix ma is complex, the Cabibbo angle is modified by 
phases, A = Aexp{z(0^ — ad)}- The resulting leptonic mixing matrix is indeed of the 
wanted form ( |130| ) with all matrix elements 0(1), except Ui 3 , 

U 13 = is 23 e^- p+ ^ - Js 23 e^ + ^ ) = 0(X, e) ~ 0.1 , (152) 
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which is close to the experimental limit. 

Let us now consider the CP violation in neutrino oscillations. Observable effects are 
controlled by the Jarlskog parameter (cf. (|84])) for which one finds 

J '-7~^ sin5 > ( 153 ) 

where 5 is some function of the unknown parameters 0(1). Due to the large neutrino 
mixing angles, J; can be much bigger than the Jarlskog parameter in the quark sector, 
J q = £>(A 6 ) ~ 1(T 5 . 

According to the seesaw mechanism neutrinos are Majorana fermions. This can be 
directly tested in neutrinoless double /3-decay. The decay amplitude is proportional to 
the complex mass 

(m) = y U^rm = W (A 2 M V 4 ^-^-^ - 2Aee^- Qd+2< «) m 3 . (154) 

1 ' ^— ' 1 + 77 r 

i 

With m 3 ~ a/ Am 2 tm ~5x 10~ 2 eV this yields (m) ~ 1CT 3 eV, more than two orders of 
magnitude below the present experimental upper bound [|14| . 



6 Leptogenesis 

One of the main successes of the standard early-universe cosmology is the prediction of 
the abundances of the light elements, D, 3 He, 4 He and 7 Li. Agreement between theory 
and observation is obtained for a certain range of the parameter 77, the ratio of baryon 
density and photon density |T0| ], 



n B 



r]= — = (1.5-6.3) x 1(T 10 , (155) 



where the present number density of photons is n 7 ~ 400/cm 3 . Since no significant 
amount of antimatter is observed in the universe, the baryon density yields directly the 
cosmological baryon asymmetry, Y# = (nB—ng)/s ~ r)/7, where s is the entropy density. 

A matter-antimatter asymmetry can be dynamically generated in an expanding uni- 
verse if the particle interactions and the cosmological evolution satisfy Sakharov's con- 
ditions [|70|1, 

• baryon number violation , 

• C and CP violation , 

• deviation from thermal equilibrium . 
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Although the baryon asymmetry is just a single number, it provides an important re- 
lationship between the standard model of cosmology, i.e. the expanding universe with 
Robertson- Walker metric, and the standard model of particle physics as well as its ex- 
tensions. 

At present there exist a number of viable scenarios for baryogenesis. They can be 
classified according to the different ways in which Sakharov's conditions are realized. 
In grand unified theories B and L are broken by the interactions of gauge bosons and 
leptoquarks. This is the basis of classical GUT baryogenesis . Analogously, the lepton 
number violating decays of heavy Majorana neutrinos lead to leptogenesis |72| . In the 
simplest version of leptogenesis the initial abundance of the heavy neutrinos is gener- 
ated by thermal processes. Alternatively, heavy neutrinos may be produced in inflaton 
decays, in the reheating process after inflation, or by brane collisions f75fl . The observed 
magnitude of the baryon asymmetry can be obtained for realistic neutrino masses |74| . 

The crucial deviation from thermal equilibrium can also be realized in several ways. 
One possibility is a sufficiently strong first-order electroweak phase transition which 
makes electroweak baryogenesis possible f75| . For the classical GUT baryogenesis and 
for leptogenesis the departure from thermal equilibrium is due to the deviation of the 
number density of the decaying heavy particles from the equilibrium number density. 
How strong this departure from equilibrium is depends on the lifetime of the decaying 
heavy particles and the cosmological evolution. 



6.1 Baryon and lepton number at high temperatures 

The theory of baryogenesis involves non-perturbative aspects of quantum field theory and 
also non-equilibrium statistical field theory, in particular the theory of phase transitions 
and kinetic theory. A crucial ingredient is the connection between baryon number and 
lepton number in the high-temperature, symmetric phase of the standard model. Due to 
the chiral nature of the weak interactions B and L are not conserved. At zero temperature 
this has no observable effect due to the smallness of the weak coupling. However, as the 
temperature approaches the critical temperature T EW of the electroweak transition, B- 



and L-violating processes come into thermal equilibrium |76 . 

The rate of these processes is related to the free energy of sphaleron-type field config- 
urations which carry topological charge. In the standard model they lead to an effective 



interaction of all left-handed quarks and leptons |52| (cf. fig. |20| ). 

O b +l = Y\_ {QLiququki) , (156) 
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Figure 20: One of the 12-fermion processes which are in thermal equilibrium in the 
high-temperature phase of the standard model. 



which violates baryon and lepton number by three units, 

AB = AL = 3 . (157) 

The evaluation of the sphaleron rate in the symmetric high-temperature phase is a 
complicated problem. A clear physical picture has been obtained in Bodeker's effective 
theory [jn| according to which low-frequency gauge field fluctuations satisfy an equation 
analogous to electric and magnetic fields in a superconductor, 

DxB=(iE-(. (158) 

Here ( represents Gaussian noise, and a is a non-abelian conductivity. The sphaleron 
rate can then be written as, 

Tsph ^ (14.0 ± 0.3)-(aJ) 5 . (159) 
a 



Lattice simulations JT8J have confirmed early estimates that B- and L- violating processes 
are in thermal equilibrium for temperatures in the range 

T EW ~ 100 GeV < T < T SPH ~ 10 12 GeV . (160) 

Sphaleron processes have a profound effect on the generation of the cosmological 
baryon asymmetry, in particular in connection with the dominant lepton number violat- 
ing interactions between lepton and Higgs fields, 

£al=2 = \h fLiWhrf + h.c. . (161) 
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Figure 21: Tree level and one-loop diagrams contributing to heavy neutrino decays. 

As discussed in section 3.3, these interactions arise from the exchange of heavy Majorana 
neutrinos, with <p — Hi and f\j = —(hvM~ l h^)ij. In the Higgs phase of the standard 



model, where the Higgs field acquires a vacuum expectation value, the interaction ( |161| ) 
leads to Majorana masses for the light neutrinos v e) and v T . 



One may be tempted to conclude from eq. ( p. 571 ) that any B + L asymmetry generated 
before the electroweak phase transition, i.e. at temperatures T > Tew, will be washed 
out. However, since only left-handed fields couple to sphalerons, a non-zero value of B+L 
can persist in the high-temperature, symmetric phase in case of a non-vanishing B — L 
asymmetry. An analysis of the chemical potentials of all particle species in the high- 
temperature phase yields a relation between the baryon asymmetry Yg = (ris — n B )/s 
and the corresponding B — L and L asymmetries Yb-l and Yl, respectively, 

Y B = a Y B -l = Y L . (162) 

a — I 

The number a depends on the other processes which are in thermal equilibrium. If these 
are all standard model interactions one has a = 8/23 in the case of two Higgs doublets 

P3- 



From eq. ( |162| ) one concludes that the cosmological baryon asymmetry requires also 
a lepton asymmetry, and therefore lepton number violation. This leads to an intriguing 
interplay between Majorana neutrinos masses, which are generated by the lepton-Higgs 
interactions ( |161|) , and the baryon asymmetry: lepton number violating interactions are 
needed in order to generate a baryon asymmetry; however, they have to be sufficiently 
weak, so that they fall out of thermal equilibrium at the right time and a generated 
asymmetry can survive until today. 



6.2 Thermal leptogenesis 

Let us now consider the simplest possibility for a departure from thermal equilibrium, 
the decay of heavy, weakly interacting particles in a thermal bath. We choose the heavy 
particle to be the lightest of the heavy Majorana neutrinos, Ni = N = N c , which can 
decay into a lepton Higgs pair l<p and also into the CP conjugate state l(f>, 

N->l<j>, N^l$. (163) 



44 



In the case of CP violating couplings a lepton asymmetry can be generated in the decays 
of the heavy neutrinos N, 

r(iv-ty) = i(i + e a )r, r(iv-Z$ = i(i- ei )r, (164) 

where T is the total decay width and £1 <C 1 measures the amount of CP violation. 




Figure 22: AL = 1 processes: decays and inverse decays of a heavy Majorana neutrino. 
The CP asymmetry E\ arises from one- loop vertex and self-energy corrections (fig. |21} 



80 , |81| , jS^] . It can be expressed in a compact form, which in the mass eigenstate basis 



of the heavy neutrinos reads, 



o Im {rrv D m v m* D ) 
e x ~ — sign(Mi) ^ — J -± . (165) 



167T ° v ±y v\m\ 
Here, in addition to the mass matrices rriD and m„, an effective neutrino mass 

m} D rri£> 



mi = V |Mi| /U , (166) 



appears, which is a sensitive parameter for successful leptogenesis [Q. Note, that the 
maximal CP asymmetry is related to the mass Mi of the heavy Majorana neutrino N\ 
|84|| . In the case of mass differences of order the decay widths, |Mj — M\\ = (^(r^ri), 
i ^ 1, the CP asymmetry is enhanced fl8~5| . Contrary to the case considered here, decays 
of A^2 or A^3 may be the origin of the baryon asymmetry [|36 |, if washout effects are 
sufficiently small. 

The generation of a baryon asymmetry is an out-of-equilibrium process which is 
generally treated by means of Boltzmann equations ||71|| . The main processes in the 
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thermal bath are the decays and the inverse decays of the heavy neutrinos (fig. |22|), and 
the lepton number conserving (AL = 0) and violating (AL = 2) processes (fig. ^). In 
addition there are other processes, in particular those involving the t-quark, which are 
important in a quantitative analysis [[87], jB3fl . 




Figure 23: AL = and AL = 2 lepton Higgs processes. 

A typical solution of the Boltzmann equations is shown in fig. Here the ratios of 
number densities and entropy density, 

Y x = — , (167) 
s 

are plotted, which remain constant for an expanding universe in thermal equilibrium. A 
heavy neutrino, which is weakly coupled to the thermal bath, falls out of thermal equilib- 
rium at temperatures T ~ M, since its decay is too slow to follow the rapidly decreasing 
equilibrium distribution ~ exp(— j3M). This leads to an excess of the number density, 
un > n e ^. CP violating partial decay widths then yield a lepton asymmetry which, by 
means of sphaleron processes, is partially transformed into a baryon asymmetry. 



The CP asymmetry e± (|165|) leads to a lepton asymmetry in the course of the cos- 



mological evolution, which is then partially transformed into a baryon asymmetry |72 
by sphaleron processes, 

n L -n T ex 

Y B = = c s k, — , (168) 

s g* 

where c$ = —8/15 is the sphaleron conversion factor in the case of two Higgs doublets. In 
order to determine the washout factor k < 1 one has to solve the Boltzmann equations. 
Compared to other scenarios of baryogenesis this leptogenesis mechanism has the advan- 
tage that, at least in principle, the resulting baryon asymmetry is entirely determined 
by neutrino properties. 

We are now ready to determine the baryon asymmetries predicted by the models of 
neutrino masses discussed in section 5. Since for the Yukawa couplings only the powers 
in e are known, we will also obtain the CP asymmetry and the corresponding baryon 
asymmetry to leading order in e, i.e. up to unknown factors 0(1). 
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Figure 24: Time evolution of the number density to entropy density ratio. AtT^M the 
system gets out of equilibrium and an asymmetry is produced. 

Consider first the model with symmetry group SU(5)xU(l)i? (cf. section 5.2). One 

(169) 



easily obtains from eqs. ( |134j ), ( 165 ) and table 

3 



5i 



16tt 



With e 2 ~ 1/300 and g* ~ 100 this yields the baryon asymmetry, 

Y b ~k 10" 8 . 



(170) 



For k ~ 0.1. ..0.01 this is indeed the correct order of magnitude. The baryogenesis 
temperature is given by the mass of the lightest of the heavy Majorana neutrinos, 



T B ~ Mi ~ e 4 M 3 ~ 10 10 GeV . 



(171) 



This is essentially the model studied in ref. [)3S}1 , where the CP asymmetry is determined 
by the mass hierarchy of light and heavy Majorana neutrinos. It is remarkable that the 
observed baryon asymmety is obtained without any fine tuning of parameters, if B — L is 
broken at the unification scale Agut- Note, that the generated baryon asymmetry does 
not depend on the flavour mixing of the light neutrinos. 

Qualitatively, the SO(10) model discussed in section 5.3 is very similar. For a class 
of parameters corresponding to the LMA MSW-solution with a, p, rj = 0(1), one finds 
for the CP asymmetry |56| , 
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(172) 



Note, that except for a, the CP asymmetry is determined by parameters of the quark 
mass matrices. Similar results have been obtained in refs. [^, |90| for different models. 
For a recent analysis within SO(10) models, which favours the Just-So and SMA solar 
neutrino solutions, see ref. 
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Figure 25: Time evolution of the heavy neutrino number density and the lepton asymmetry 
for the SO (10) model. The upper solid line shows the solution of the Boltzmann equations 
for the right-handed neutrinos; the equilibrium distribution is represented by the dashed 
line. The absolute value of the lepton asymmetry Yl is given by the lower solid line, 
and the hatched area shows the lepton asymmetry corresponding to the observed baryon 



asymmetry 



O 




Figure 26: Same as fig. E3I with initial asymmetry Y[ n ~ 10 [92| 
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Numerically, with e ~ 0.1 one finds e\ ~ 10~ 7 and 

\M X \ ~ e e (1 + |r/|2K2 ~ 10 9 GeV , (173) 
crm 3 

with rhi ~ {\v\ 2 + IpI 2 )^/^! + M 2 )) ~ 10~ 2 eV. The baryon asymmetry is then given 

by 

Y B ~ —K sign(o-) sin (<p u - a u ) x 1(T 9 . (174) 

The solution of the full Boltzmann equations is shown in fig. The initial condi- 
tion at a temperature T ~ 10Mx is chosen to be a state without heavy neutrinos and 
with vanishing lepton asymmetry. The Yukawa interactions are sufficient to bring the 
heavy neutrinos into thermal equilibrium. At temperatures T ~ M\ the familiar out-of- 
equilibrium decays set in, which leads to a non-vanishing baryon asymmetry. The dip 
in fig. ^5] is due to a change of sign in the lepton asymmetry at T ~ M\. The final 
asymmetry is about 1/3 of the observed value, which lies within the present range of 
theoretical uncertainties. 

A very important question for leptogenesis, and baryogenesis in general, is the de- 
pendence on initial conditions. As fig. [25| demonstrates, the heavy neutrinos are initially 
indeed in thermal equilibrium. One may also wonder, how sensitive the final lepton 
asymmetry is to an initial asymmetry which may have been generated by some other 
mechanism. The SO(10) model under consideration turns out to be very efficient in 
establishing a symmetric initial state. This can be seen in fig. |2T^ where the maximal 
asymmetry Y[ n ~ 1CT 2 has been assumed as initial condition. Within one order of mag- 
nitude in temperature this initial asymmetry is washed out by eight orders of magnitude! 
Hence, the final baryon asymmetry is a definite prediction of the theory, independent of 
initial conditions. 

In summary, the experimental evidence for small neutrino masses and large neutrino 
mixings, together with the known small quark mixings, have important implications for 
the structure of grand unified theories. In SU(5) models this difference between the lepton 
and quark sectors can be explained by U{1)f family symmetries. In these models the 
heavy Majorana neutrino masses are not constrained by low energy physics, i.e. light 
neutrino masses and mixings. Successful leptogenesis is possible, but it depends on the 
choice of the heavy Majorana neutrino masses. 

In SO (10) models the implications of large neutrino mixings are much more stringent 
because of the connection between Dirac neutrino and up-quark mass matrices. The 
requirement of large neutrino mixings then determines the relative magnitude of the 
heavy Majorana neutrino masses in terms of the known quark mass hierarchy. This leads 
to predictions for neutrino mixings and masses, CP violation in neutrino oscillations 
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and neutrinoless double /3-decay. It is remarkable that the predicted order of magnitude 
of the baryon asymmetry is also in accord with observation. 

I would like to thank the participants of the school for stimulating questions and the 
organizers for arranging an enjoyable and fruitful meeting at Beatenberg. I have also 
benefited from discussions with P. Di Bari, M. Pliimacher, A. Ringwald, S. Stodolsky 
and D. Wyler, and I thank S. Wiesenfeldt for providing some of the figures. 
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